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SUMMARY 

The Bozzolo-Ferrante-Smith (BFS) method for alloy energetics is applied to the study of ternary 
additions to NiAl. A description of the method and its application to alloy design is given. Two 
different approaches are used in the analysis of the effect of Ti additions to NiAl. First, a thorough 
analytical study is performed, where the energy of formation, lattice parameter and bulk modulus 
are calculated for a large number of possible atomic distributions of Ni, A1 and Ti. Substitutional 
site preference schemes and formation of precipitates are thus predicted and analyzed. The second 
approach used consists of the determination of temperature effects on the final results, as obtained 
by performing a number of large scale numerical simulations using the Monte Carlo-Metropolis 
procedure and BFS for the calculation of the energy at every step in the simulation. The results 
indicate a sharp preference of Ti for A1 sites in Ni-rich NiAl alloys and the formation of ternary 
Heusler precipitates beyond the predicted solubility limit of 5 at. % Ti. Experimental analysis of 
three Ni-Al-Ti alloys confirms the theoretical predictions. 

KEYWORDS: Alloys, Intennetallics, BFS Method, Defect Structure, Semiempirical Methods, 
Ternary Alloy, Heusler phase 
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INTRODUCTION 


The development of more efficient gas turbine engines and power plants for future supersonic 
transports like the high speed civil transport (HSCT) depends upon the advancement of new high- 
temperature materials with temperature capabilities exceeding those of current nickel-base super- 
alloys. Ordered intermetallic alloys, in particular the B2 structured NiAl, have long been consid- 
ered a prime candidate to replace superalloys in the combustion and turbine sections of aircraft 
engines due to a number of property advantages including the recent development of alloying 
schemes for enhanced creep strength [1], 

One of the most basic of these strategies is to combine single crystal processing and alloying 
with reactive elements such as Ti, Hf, Zr and Ta in order to generate extremely creep resistant 
materials [2]. This approach has been so effective that single-crystal NiAl alloy turbine vanes 
have been successfully engine tested [3], Not only did the parts survive engine testing but they 
displayed superior performance to superalloy components undergoing the same evaluation. How- 
ever, alloy design, while successful to this point, has been entirely empirical in nature. Further- 
more, limited understanding of the basic microstractures and structure-property relationships has 
occurred due to this empirical alloy design approach. For example, Ti additions on the order of 
2. 5-3.0 at. % result in a 200-5000 fold reduction in creep rate over binary NiAl single crystals and 
Hf additions can have an even larger effect, but the basic strengthening mechanisms and even the 
equilibrium microstructures of these materials are unclear [4], 

While it is well known that ternary and higher order alloying additions to structural interme- 
tallics are essential for the optimization of physical, chemical, structural, and mechanical proper- 
ties, it is only recently that substantial theoretical work has been directed towards a thorough 
understanding of the atomic processes involved. In the simple cases of FeAl and NiAl, there are 
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abundant, sometimes conflicting results from a variety of techniques, providing some guidance as 
to the defect structures in these binary compounds as a function of stoichiometry [1,5-7]. From a 
theoretical standpoint, first-principles calculations have also provided very valuable information 
regarding this issue [8]. However, a full understanding of the atomic-level processes involved 
and a complete description of the composition and temperature dependence of point defects and 
site substitution schemes for alloying additions are not yet available for NiAl and FeAl. This 
problem is even more pronounced for alloy systems which have not been as thoroughly studied 
and is nearly intractable for quaternary and higher order additions to an ordered A-B system. 

While it would be extremely advantageous to perform accurate quantitative and qualitative 
theoretical analysis of these multicomponent systems, this is beyond the scope of most currently 
available techniques. Theoretical predictions of alloy properties are generally obtained by two 
approaches, first-principles and semiempirical. First-principles methods involve solving Schro- 
edinger’s equation for the system at hand, while semiempirical methods attempt to develop 
approaches to model energetics with empirical input to obtain unknown parameters. Ideally, first 
principles calculations are best suited for providing the most accurate and consistent framework 
for such studies. However, the small differences in energy related to basic issues like preferential 
site substitution impose the need for intensive electronic structure calculations. Moreover, the 
complexity of the multi-component systems imposes further limitations, and thus, not enough 
information can be realistically expected from a first-principles approach because the CPU inten- 
sive calculations must be performed for many possible, large configurations. This limits the use- 
fulness of first principles approaches as economical predictive tools. 

On the other hand, semiempirical approaches partially solve this problem by introducing a 
much needed degree of efficiency, which enables one to look at larger systems, at the expense of 
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a hopefully minimal loss of physical accuracy. Most importantly, they are useful in providing a 
global view of the mechanisms or processes at hand. However, their range is limited, in most 
cases, to a few elements or some particular crystallographic structure. Furthermore, methods that 
exhibit a great deal of accuracy for bulk materials fail to reproduce some of the most basic proper- 
ties of surfaces or grain boundaries, hi addition, their field of applicability to monatomic systems 
or very specific binary systems (for which specific parameterizations or potentials are often devel- 
oped) is limited due to the lack of transferrability of the parameters used. Recently, a new 
semiempirical technique was developed which avoids these current limitations. The BFS (Boz- 
zolo- Ferrante-Smith) method for alloys [10] is particularly designed to deal with complex sys- 
tems and geometries; it has no constraints regarding the atomic species under consideration, their 
number, or the crystallographic structure. The BFS method provides a simple framework for large 
scale computer simulations with the appropriate formal background, but without the limitations 
that exist in competing techniques. 

After a brief introduction of the BFS method, this paper concentrates on its application to the 
study of ternary Ni-Al-Ti alloys, focusing on physical properties, resulting defect structures, the 
solubility of Ti in NiAl and precipitation of second phases. A discussion of the lattice parameter 
dependence with composition based on a simple, approximate expression is also included. More- 
over, the results include Monte Carlo temperature-dependent large- scale simulations providing 
some insight on final equilibrium microstructures. Finally, owing to the engineering significance 
of NiAl-Ti alloys and to develop confidence in the modeling effort, the theoretical results derived 
from the BFS method are compared to the results from a concurrent transmission electron micros- 
copy study. 
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THE BFS METHOD 


Since its inception a few years ago, the BFS method [9-11] has been applied to a variety of 
problems, ranging from bulk properties of solid solutions of fee and bcc binary alloys [11-12] to 
more specific applications like the energetics of bimetallic tip- sample interactions in an atomic 
force microscope [13]. More recently, the BFS method was used to deal with alloy surfaces as 
well as the phase stability in bulk alloys, which has been illustrated in detailed studies and simula- 
tions of surface alloys [14] and the design of Ni-based quaternary alloys [10], providing a founda- 
tion for the work presented in this paper. In addition, through BFS we can derive simple, 
approximate expressions that describe trends in segregation as well as elucidating driving mecha- 
nisms for these phenomena [15]. Also, simple expressions can be derived for predicting the com- 
position dependence of bulk alloy properties based solely on pure component properties (the BF 
rule) [16]. 

In what follows, we provide a brief description of the operational equations of BFS in order to 
introduce some concepts used later on. The reader is encouraged to seek further details in previ- 
ous papers where a detailed presentation of the foundation of the method, its basis in perturbation 
theory, and a discussion of the approximations made are clearly shown [9,10] . 

The BFS method provides a simple algorithm for the calculation of the energy of formation of 
an arbitrary alloy (the difference between the energy of the alloy and that of its individual constit- 
uents). hi BFS, the energy of formation is written as the superposition of individual contributions 
of all the atoms in the alloy 

AH = £e, ( 1) 
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For each atom, we partition the energy into two parts: a strain energy and a chemical energy 
contribution. The strain energy is based on the atomic positions of the neighboring atoms to atom 
i, regardless of their chemical identity. For its calculation, we use the actual geometrical distribu- 
tion of the atoms surrounding atom /, computed as if all of its neighbors were of the same species 
as atom i. In this sense, the BFS strain energy differs from the commonly defined strain energy in 
that the actual chemical environment is replaced by that of a monatomic crystal. Its calculation is 
then straightforward, even amenable to first-principles techniques. In our work, we use Equivalent 
Crystal Theory (ECT) [17] for its computation, due to its proven ability to provide accurate and 
computationally economical answers to most general situations. In all cases considered in this 
work, a rigorous application of ECT is reduced to that of its two leading terms, describing average 
density contributions and bond-compression anisotropies and neglecting the three- and four-body 
terms dealing with the bond angle and face-diagonal anisotropies [17]. 

The chemical environment of atom / is considered in the computation of the BFS chemical 
energy contribution, where the surrounding atoms are forced to occupy equilibrium lattice sites 
corresponding to the reference atom ;. Building on the concepts of ECT, a straightforward 
approach for the calculation of the chemical energy is defined, properly parameterizing the inter- 
action between dissimilar atoms. 

Thus defined, the BFS strain and chemical energy contributions take into account different 
effects, i.e. geometry and composition, computing them as isolated effects. A coupling function g 
restores the relationship between the two terms, defined in such away as to properly consider the 
asymptotic behavior where chemical effects are negligible for large separations between dissimi- 
lar atoms. 
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Simply, the contribution from each atom i to the energy of formation A H is 


8,. = e s + g e c = e s + g(e c - e c ) 


( 2 ) 


where e c is the chemical energy of the atom (e c ) referenced to its pure ground state crystal (e c ). 
The strain energy contribution is obtained by solving the ECT perturbation equation 


l 

NR[e~ aRl + MR p 2 e ‘ £ 


p -(a + S(r ))r ■ 

r j e 


(3) 


where N and M are the number of nearest- and next-nearest neighbors respectively, and where p, /, 
a and X are ECT parameters that describe element / (see Ref. 17 for definitions and details), r 
denotes the distance between the reference atom and each of its neighbors and S(r) describes a 
screening function [17] and the sum runs over nearest and next nearest neighbors. This equation 
determines the lattice parameter a s of a perfect, equivalent crystal where the reference atom i has 
the same energy as it has in the geometrical environment of the alloy. R / and R 2 denote the near- 
est- and next-nearest neighbor distances in this equivalent crystal. 

Once the lattice parameter of the (strain) equivalent crystal , a s , is determined, the BFS strain 
energy contribution is computed using the Universal Binding Energy Relation (UBER) of Rose et 
al. [18], which contains all the relevant information concerning a single-component system: 

(4) 


e s = E c 1 - (1 + a s )e 


where E c is the cohesive energy of atom / and where the scaled lattice parameter a* is given by 


* (« s ~ a e) 

a S = ( 1 ] 


( 5 ) 
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where q is the ratio between the equilibrium Wigner-Seitz radius and the equilibrium lattice 
parameter a e . 

The BFS chemical energy is obtained by a similar procedure. As opposed to the strain energy 
term, the surrounding atoms retain their chemical identity, but are forced to be in equilibrium lat- 
tice sites of an equilibrium (otherwise monatomic) crystal i . The BFS equation for the chemical 
energy is given by 


NR\'e 


v/R , 


+ MR 2 ‘e 


= x 


"i/i* 


+ M i/i e 




(6) 


where N jk and M jk are the number of nearest- and next-nearest neighbors of species k to atom i. 
The chemical environment surrounding atom i is reflected in the parameters a jk , given by 


a ik = a i + A ki 


(7) 


where the BFS parameters A ki (a perturbation on the single-element ECT parameter o.J describe 
the changes of the wave function in the overlap region between atom i and k. Once Eq. (6) is 
solved for the equivalent chemical lattice parameter a c , the BFS chemical energy is then 


e 


c 



1 - (1 + ci c )e 



( 8 ) 


where Y; = 1 if a c > o and y,- = - 1 if a* c < o . The scaled chemical lattice parameter is given by 
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A reference chemical energy contribution e c is also computed, in order to free the BFS chem- 
ical energy from any structural information [9,10]. This is accomplished by repeating the process 
for computing e c (using Eqs. 6 and 8) but setting the BFS parameters A ki to zero. The total BFS 
chemical energy e c (Eq. 2) is then 

e c = e c - e c 0 • ( 10) 

Finally, as mentioned above, the BFS chemical and strain contributions are linked by a cou- 
pling function g which describes the influence of the geometrical distribution of the surrounding 
atoms in the relevance of the chemical effects, and is given by 



(11) 


In this work we used the BFS interaction parameters A determined following the procedure 
outlined in Ref. 9. These parameters are obtained from first-principles, all-electron, density func- 
tional calculations of the elemental constituents and ordered binary compounds of these elements 
[19,20]. The particular implementation used is the Linear- Muffin-Tin Orbitals (LMTO) method 
[10] in the Atomic Sphere Approximation. As mentioned above, in order to provide parameters to 
the BFS method, we need to calculate the equilibrium properties of the elemental solids in the 


NAS A/TM— 1999-208892 


9 



same symmetry as the compound to be studied, since BFS is referenced to the ground state prop- 


erties of the system in that symmetry. Therefore, for this study, Ni, A1 and Ti input parameters 
were calculated for the bcc symmetry (Table 1). Also, BFS interaction parameters A ab and A ba 
were determined for Ni-Al, Ni-Ti and Ti-Al in the bcc symmetry by fitting to LMTO-computed 


TABLE 1 

COMPUTED CONSTANTS AND INPUT FOR ECT 



LMTO results 



ECT Parameters 

Atom 

Lattice 

Parameter 

O 

(A) 

Cohesive 

Energy 

(eV/atom) 

Bulk 

Modulus 

(GPa) 

P 

a 

(A' 1 ) 

1 

O 

(A) 

X 

O 

(A) 

Ni 

2.752 

5.869 

249.2 

6 

3.0670 

0.763 

0.2716 

A1 

3.192 

3.942 

77.3 

4 

1.8756 

1.038 

0.3695 

Ti 

3.213 

6.270 

121.0 

6 

2.6805 

1.048 

0.3728 


LMTO results for the lattice parameter, cohesive energy and bulk modulus for the bcc phases of Ni, A1 and 
Ti. The last four columns list the ensuing ECT parameters determined from the LMTO results. 


TABLE 2 
INPUT FOR BFS 



o 

BFS Parameters (A 

- 1 ) 

A-B 

Aab 

Aba 

Ni-Al 

-0.05813 

0.0822 

Ni-Ti 

-0.06587 

0.4610 

Al-Ti 

-0.06360 

0.2283 


BFS parameters A ab and A ba for Ni-Al, Ni-Ti and Al-Ti determined by fitting the lattice parameter and 
energy of formation of the corresponding B2 compounds via LMTO calculations. 
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energies of formation AH and equilibrium lattice parameter a 0 for the corresponding B2 alloys 
(Table 2). 

Once these parameters are computed, they remain the same for any other calculation involving 
any of these elements, requiring no further adjustment or replacement unless the crystal symmetry 
we are dealing with changes. 

Ni-Al-Ti ALLOYS 

From a practical standpoint two-phase alloys based on a NiAl matrix reinforced by Heusler 
(Ni 2 AlTi) precipitates have been receiving increased experimental attention due to their potential 
as high temperature structural alloys [2-4,25-28]. The Heusler phase has a cubic L2 1 structure, in 
which the unit cell is comprised of eight bcc unit cells with A1 and Ti atoms occupying two sets of 
octahedral sites located at body-center locations (Fig. 1). It is related to the B2-NiAl structure in 
that every other A1 site in the NiAl lattice becomes occupied by a Ti atom in an ordered fashion. 



a 2 bc 


• a Hb Oc 


Fig. 1 : Schematic illustration of the Heusler phase and its relation to B2 compounds. 


NAS A/TM— 1999-208892 


11 


Because of the similarity in lattice structure between B2 and L2|, it is possible to develop precip- 
itation hardened alloys similar to y/y’ nickel base superalloys. However, in order to design within 
this new family of B2/L2 1 alloys (also known as (3/(3’ superalloys), microstructural factors such as 
the solubility of Ti in NiAl, the defect structure within the solid solution NiAl-Ti phase, the lattice 
misfit between the two phases, and eventually the effect of quaternary and higher order alloying 
additions on all of these factors need to be determined in detail. However, even with the amount of 
experimental work performed on the ternary NiAl-Ti system over the last decade [2-4,25-28], 
details of these microstructural features are still vague while those for quaternary and higher order 
systems are completely lacking. 

Fu et al. [29] have recently presented a comprehensive first principles study of site substitution 
on both the defect structure of FeAl and NiAl compounds as well as the behavior of ternary (Ti or 
Cr in FeAl and Fe in NiAl) additions to these systems. Unfortunately, no NiAl+Ti results were 
included. In a much earlier study, Tso and Sanchez [30] performed some thermodynamic model- 
ing of the Ni-Al-Ti system but the results focused on modeling of the various binary phase dia- 
grams with limited attention and no detail presented concerning the case of ternary alloys. 

In this paper, we concentrate on Ni-rich NiAl alloys with Ti alloying additions, prevalent in 
industrial applications [3,4], It is essential that the basic properties of the NiAl alloy be accurately 
reproduced by the use of the BFS method. As seen in Ref. 9, the defect structure of binary NiAl 
alloys consists primarily of substitutional antisite defects on the Al sublattice and as a result the 
lattice parameter monotonically decreases with increasing Ni concentration [1], In a previous 
application of the BFS method to the study of the zero temperature defect structure of NiAl alloys 
[9,31], the computational results correctly identified the energetically favored defect structures, 
including the correct substitutional defect scheme in Ni-rich alloys. It was also shown that while 
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Fig. 2 : Lattice parameter of non-stoichiometric NiAl alloys as a function of Ni concentration, normal- 
ized to their stoichiometric values. The solid circles denote experimental results from different investi- 
gators (see Ref. 1). The lines denote the BFS predictions. 


Ni vacancies are mostly responsible for the observed behavior in Al-rich alloys, evidence was 
found for the possibility of vacant sites in both sublattices as well as a trend for clustering of 
vacancies. Both features have been observed in recent experiments [32], 

In particular, the quantitative accuracy of the BFS results for binary NiAl is highlighted by the 
predicted values of the lattice parameter. A survey of available experimental data indicates that the 
lattice parameter in the Ni-rich region varies linearly with concentration [1], Correspondingly, the 
BFS results can be adjusted to a similar linear regime with a very small departure from the exper- 
imental expressions: 


a_ \ 

a 0 'Exp 


1.0384- 0.00079 lAx m 


( 12 ) 
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and 


(— ) = 1.0396 - 0. 0007689a w (13) 

^ a O'BFS 

where x Nj denotes the atomic concentration of Ni. Fig. 2 displays both the experimental values as 
well as the two linear regions representing the theoretical results. The qualitative and quantitative 
agreement with experiment certify the validity of the Ni-Ni, Al-Al and Ni-Al BFS parameters 
used in this work (listed in Table 2). From a methodological point of view, the good agreement 
between theoretical and experimental results that follow can be taken as evidence for the validity 
of the BFS parameters used in this paper. 

Analytical Structures 

By exploiting the computational simplicity of semiempirical methods and in particular the 
BFS method, it is possible to analize the energetics of a large number of Ni-Al-Ti alloys. These 
alloys represent a wide range of concentrations and different atomic distributions with different 
degrees and types of ordering. For reasons of simplicity and to enable comparison with available 
experimental data on industrially relevant compositions, we restrict our calculations to Ni-rich 
Ni 50 Algo_ x Ti x alloys. The large size of the necessary computational cells combined with the 
number of different elements considered results in a large number of possible distributions of 
these atoms within the computational cell chosen. Consequently, an efficient and economical - but 
physically sound - technique is required to examine all pertinent configurations of atoms. While a 
search of every possible atomic arrangement in a reasonably large computational cell is not 
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impossible, we seek to develop a simplified approach that would not only make it more physically 
intuitive, but also amenable to extensions to larger systems without fearing reaching the limit of 
the available computational tools. 

In this section, we use a similar approach to that used in Ref. 9, representing a compromise 
between a detailed analytical study and an exhaustive statistical survey for each concentration 
and atomic distribution. The procedure followed in Ref. 9 consisted of definining a specific set of 
atomic 'configurations’, and then using the BFS method to compute the energy of formation of 
each configuration, as well as its equilibrium lattice parameter and bulk modulus. A large set of 
high symmetry configurations covering a wide composition range was defined in this maimer. 
The contents of this catalogue of configurations is, to great extent, arbitrary. However, a ternary 
system is still simple enough to develop an intuitive algorithm for constructing configurations of 
potential interest. A similar approach is used hi this work to examine ternary alloys, building on 
the originally defined set of configurations including some selected arrangements of atoms that 
represent different ordering patterns. The purpose of this approach is to 1) relate general trends of 
the lattice parameter and the energy of formation to changes in concentration and atomic distribu- 
tion and 2) identify metastable structures, that is, configurations with energy close to that of the 
ground state, which might have a high probability of appearing in the actual alloy depending on 
the conditions prevalent during its processing. Therefore, this 'catalogue’ of alternative configura- 
tions serves the purpose of identifying ordering trends and provides physical insight to the micro- 
structure of the actual alloy. Since the number of possible configurations could become 
prohibitively large, we believe that this selection procedure, when properly undertaken, is suffi- 
cient to gam an adequate understanding of the alloy composition under investigation. Some 
selected configurations are described in Fig. 3. The labeled 72-atom computational cell is shown 
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Fig. 3 (a) The computational cell used to generate the catalogue shown in the Appendix, (b-c) Two exam- 
ples of the configurations of NiAl-Ti alloys used in this calculation, corresponding to (b) x Ti = 13.89 and 
(c) x Ti = 25.00. In these last two cases, the label indicates the type of substitution. X n -> Y m indicates an 
atom of species X, originally in site n, substituting for an atom of species Y in site m. 
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ill Fig. 3. a to clarify the atomic exchanges that define the different configurations in the set. Two 
examples of such configurations are shown in Fig. 3.b and 3.c, and a complete list of the configu- 
rations studied in this work is included in the Appendix. 

The obvious difficulty in defining an appropriate set of configurations for the study at hand is 
the lack of a priori knowledge of the ground state structure for a given composition or tempera- 
ture. In some cases (i.e, low concentration solutes) the ground state might have easily predictable 
features but, in general, no guarantee exists that the ground state configurations for each concen- 
tration are included. However, a close examination of the results usually give a clear indication 
whether a particular configuration could be 'missing’ from the set if, for example, there are data 
points ‘missing’ in what would otherwise be an easily recognizable trend. 

In this work, the results for a set of over 150 different configurations or arrangements of a 72- 
atom computational cell are shown and discussed. Fig. 3 shows two of these configurations and 
introduces the notation used in the Appendix to identify the positions of Ni, A1 and Ti atoms in the 
lattice. Obviously, if this set was complete and if it included the ground state for a certain compo- 
sition, the energy of formation for such a configuration would be lower than any other with the 
same concentration of elements. Moreover, if the catalogue includes the ground state configura- 
tions for all concentrations, it would be easy to determine the symmetries that characterize such a 
state. 

Fig. 4 summarizes the results of this type of analysis for the 150+ set of configurations: the 
energy of formation at zero temperature for Ni 5 QAl 5 o_ x Ti x alloys, for x between 0 (B2 NiAl) and 
50 at. % Ti (B2 NiTi), is shown. A close examination of the atomic distributions considered for 
each concentration -detailed in the Appendix- reveals that most of the likely ordering patterns for 
a given composition are included in this survey. Moreover, Fig. 4 shows, for each concentration, 
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'clusters’ of states very close in energy that correspond, in terms of the atomic distributions, to 
related ordering patterns. This behavior is clear throughout the whole range of compositions 
included, where both end points (NiAl and NiTi) form B2 compounds. 

If we focus on those configurations with 0< x Ti < 25 %, it is clear that a specific group of con- 
figurations, corresponding to the lowest energy states for each concentration above x Ti = 5%, fol- 

(a) 

-0.35 

-0.40 

-0.45 

a 

o 
te 

> -0.50 

<D 

X -0.55 

< 

-0.60 
-0.65 

0 10 20 30 40 50 

x Ti 

Fig. 4. Energy of fonnation (in eV/atom) of the cells listed in the Appendix for 0 < x Ti < 50. The circles 
denote those configurations characterized by Heusler ordering, while the solid squares include a variety of 
short-range order patterns as well as disordered states. 
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lows a specific ordering trend. We will return to this issue later, once some basic features 


regarding site preference schemes for the alloying addition are determined. 


Site preference ofTi in NiAl alloys 


In general, the results in Fig. 4 indicate that those configurations where Ti atoms occupy sites 
on the A1 sublattice are consistently lower in energy. The site preference of Ti atoms can be ana- 


TABLE 3 

SITE PREFERENCE IN A Ni 50 Al 50 _ x Ti x ALLOY 


x Ti 

at. %Ti 

Description of the configuration 

Energy of 
Formation 
(eV) 

Lattice 

Parameter 

(A) 


Ti(Ni)Al NN 

-0.41320 

2.870 

1.39 

Ti(Ni)Al far 

-0.38052 

2.870 


Ti(Al) 

-0.60839 

2.853 


2Ti(Ni)NN+2Ni(Al) NN 

-0.32607 

2.881 


2[Ti(Ni)Al] NN 

-0.32873 

2.882 


Ti(Ni)Al NN +Ti(Al) 

-0.36893 

2.879 

2.78 

[Ti(Ni)Al+Ti(Al)] far(2a) 

-0.37773 

2.878 

2Ti(Al) NN 

-0.59918 

2.859 


2Ti( Al) far( jYa ) 

-0.60350 

2.858 


2Ti(Al) Heusler 

-0.60350 

2.858 


2Ti( Al) far( JYa ) 

-0.60363 

2.858 


Energy of formation (in eV/atom) and lattice parameter (in A) for specific atomic configurations of Ti 
atoms within a NiAl lattice at two different concentration levels. 
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lytically determined with the BFS method by considering just a few selected configurations 
where, for a fixed concentration, Ti atoms are located in specific sites. Figure 5 shows the config- 
urations considered in such a search and Table 3 lists the corresponding values for the energy of 
formation and equilibrium lattice parameter. 

The notation used to indicate site preference is self-explanatory and is convenient for repre- 
senting more complicated situations as will be seen in subsequent studies dealing with 4- and 5- 
element systems. The notation A(B) indicates a substitutional defect where atom A occupies a B 
site. The designation A(B)C indicates an A atom occupying a B site with the displaced B atom 
occupying a C site. When necessary, the physical separation between two point defects is noted : 
for example, A(B)C NN indicates that the A atom in the B site and the B atom in the C site are 
nearest neighbors. In Table 3, the subindex far(d) denotes that the corresponding atoms are at a 
distance d greater than second neighbor distances, while the subindex Heusler indicates that the 
atoms locate themselves in A1 sites following the Heusler structure (atoms in any pair of Heusler 
sites are located in opposite comers of a face of the elementary cube in the A1 sublattice). 

The cases considered for the study of site substitution for low Ti concentration include two 
different values: 1.39 and 2.78 at. % Ti. For the first case, which represents one Ti atom in the 72- 
atom cell, two basic configurations are possible: the Ti atom occupying an A1 site (Ti(Al)) (fig. 
5. a), and the Ti atom occupying a Ni site with the displaced Ni atom located in an A1 site 
(Ti(Ni)Al). In this second type of configuration, the antistructure Ni atom can be a nearest neigh- 
bor of the Ti atom (Ti(Ni)Al^, fig. 5.b) or not ( Ti(Ni)Al^ /n fig. 5.c). The corresponding values 
for the energy of formation are listed in Table 3, clearly showing the preference of a Ti atom for 
an A1 site, A similar situation is observed at a higher Ti concentration (figs. 5.d-5.h), which corre- 
sponds to two Ti atoms in the 72-atom cell. Several possibilities now exist concerning the relative 
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(a) Ti(Al) 


(b) Ti(Ni)Al NN 





(g) 2Ti(Al)f ar 


(h) 2Ti(Al) Heusler 


Ni 



Ti 


Fig. 5 Configurations used for the calculation of Ti site preference in NiAl alloys, (a) A Ti atom in an At 
site, (b) a Ti atom in a Ni site, with the displaced Ni atom occupying a site in the At sublattice at nearest- 
neighbor distance from the Ti atom, (c) same as before but with the Ni atom located at a different site in the 
lattice, at a distance greater than next-nearest-neighbor distance, etc.. In all cases, A(B) indicates an atom 
A in a B site while A(B)C indicates an atom A in a B site with the displaced B atom in a C site. 
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location of the additional Ti atoms: they could both be located in A1 sites (2Ti(Al)), or create an 
antistructure Ni pair by occupying Ni sites with the Ni atoms occupying A1 sites (2Ti(Ni)Al), or a 
combination of these defects with several options for the relative location of the antistructure and 
substitutional atoms with regard to each other ( [Ti(Al) + Ti(Ni)Al] d ). 

All these possibilities are included in Fig. 5 and the corresponding formation energies are 
listed in Table 3. Once again there is a clear preference for Ti to occupy A1 sites compared to any 
other substitution scheme. Moreover, the small difference in energy of formation between the case 
when the two Ti atoms are nearest neighbors and the case where they are truly separated can be 
taken as the reason why at low Ti concentrations there is a close competition between a solid solu- 
tion alloy and an ordered arrangement of Ti atoms on the Al-sublattice, leading to the formation 
of a precipitate phase, as seen in Fig. 4.b (cases involving ordering of the Ti atoms are indicated 
with circles while more random configurations of the Ti atoms are indicated with solid squares). 
The preference of Ti for A1 sites is observed for the whole range of Ti concentrations considered 
in this work, leading eventually to the formation of the Heusler (Ni 2 AlTi) alloy, where Ti atoms 
would be exclusively located in the A1 sublattice of a B2 compound, leading to the L2 j structure. 

Figure 6 is an alternative way of displaying the results of the site preference calculations, in 
the fonn of an energy spectrum. The advantages of visualizing the results in this way will become 
apparent when dealing with 4- and 5-element alloys, clearly indicating the interaction between the 
different alloying additions. The first column shows the energy 'levels’ corresponding to the dif- 
ferent substitution schemes for just one Ti atom in a 72 atom cell. The second and third columns 
show results for two (x Ti =2.78) and five (x Ti =6.94) atoms respectively. The splitting of the Ti(Al) 
level for x Ti =2.78 is due to the different relative locations of the two Ti atoms: they could be non- 
interacting (' far’), situated in a Heusler-like pattern (opposite comers of a cube in the A1 sublat- 
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Energy of Formation (eV/atom) 


1.39 at % Ti 


2.78 at. % Ti 


6.94 at. % Ti 


-0.30 -l 


2Ti(Ni)Al fai Other xTi(Ni)Al+yTi(Al) 

E form > -0.30 eV/at 

2Ti(Ni)Al NN 


-0.40 


Ti(Ni)Al+Ti(Al) 

Ti(Ni)Al NN 

far 


Ti(Ni)Al+Ti(Al) 

Ti(Ni)Al NN far 


4Ti(Al)+Ti(Ni) 


-0.50 -A 


-0.60 


Ti(Al) 


2Ti(Al) 

NN 


2Ti(Al) Heusler 

2Ti(Al) 

rar 


5Ti(Al) 

far 


5Ti(Al) Heusler 


Fig. 6: Energies of fonnation (in eV/atom) for 72 atom cells containing 36 Ni atoms, 36-n A1 atoms and n 
Ti atoms (n = 1, 2, 5). The different energy states correspond to different substitutional and defect schemes, 
as indicated in Fig. 5. 
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tice, indicated with circles in Fig. 4.b) or in such a way that the two Ti atoms in the A1 sublattice 
are at next-nearest-neighbor distance from each other. The energy difference between these two 
states is extremelly small at this concentration but it still favors a random distribution of Ti atoms 
in the A1 sublattice. For x Ti =6.94 the splitting in energy levels is much more pronounced since 
many more possibilities exist for the placement of the Ti atoms. For clarity, we only show two 
states in the ground state region: the one corresponding to the Ti atoms in solid solution ('far 7 
from each other in the A1 sublattice) and the Ti atoms following a strict Heusler ordering pattern. 
The reversal in energy levels between x Ti =2.78 and x Ti =6.94 for the solid solution and the Heu- 
sler ordering options clearly indicates that the formation of Heusler precipitates is favored at 
higher Ti concentrations. Finally, the stability of the Heusler phase (Fig. 4), where Ti atoms are 
located exclusively in A1 sites, indicates that the site preference observed at low concentrations is 
the same for the whole range of concentrations for which this phase exists, consistent with the fact 
that the energy gap between Ti(Al) and Ti(Ni)Al substitutions remains basically unchanged with 
increasing Ti concentration. 

Evolution of second-phase structures within NiAI-Ti alloys 

Figure 3 introduces some selected configurations from the complete set used in this work, as 
well as the notation used hi labeling them. A complete list of configurations is included in the 
Appendix and the energy of formation computed via BFS for all the configurations is shown in 
Fig. 4. As mentioned before, one feature in Fig. 4 becomes immediately apparent: beyond 5 at. % 
Ti, there is a clear separation in formation energy between a selected type of configurations 
(denoted with circles in Fig. 4) and the rest (denoted with solid squares). The selected group of 
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configurations correspond to a particular type of ordering, where Ti atoms locate themselves 
exclusively in A1 sites in such a way that they always have A1 atoms as next nearest neighbors. 
This ordering pattern, when extended to the case when the concentration of Ti is 25 %, corre- 
sponds to the L2j ordered structure or Heusler phase. At that particular concentration, the energy 
gap between the Heusler structure (solid circle in Fig. 4.b) and any other configuration is the larg- 
est, clearly indicating the stability of the L2j phase at a stoichiometry equivalent to Ni 2 AlTi. 

Below 5 at.% Ti, those configurations representing Ti in solid solution within the matrix are 
energetically favored, however so slightly, over those where short range order dominates. This is 
more clearly demonstrated in Figure 6. This situation is reversed at a composition near 5 at. % Ti, 
with an ever more distinguishable preference for Heusler ordering against any other option with 
increasing Ti content. In other words, Heusler-like ordering becomes clearly preferred beyond a 
certain critical value somewhere near 5 at. % Ti, in spite of the fact that even at lower concentra- 
tions those configurations with short-range Heusler-type ordering are also very low in energy. 
This crude way of determining the solubility limit of Ti in NiAl, which based on these results 
could be set somewhere near 5 at. % Ti, not only establishes such a critical value but it also pro- 
vides some insight on the behavior of the system for a wide range of concentrations surrounding 
the solubility limit. 

One clear reason for the stability of the Heusler phase is the type of nearest neighbor bonds 
present. The Heusler ordering maximizes the number of energetically favorable Ni-Al and Ni-Ti 
bonds - both B2 compounds - as seen in Table 4. In addition, the second neighbor bonds are 
either Ti-Al, also energetically favorable, or Ni-Ni. The Ni-Ni bonds, due to the closeness 
between the lattice parameter of the Heusler phase to the equilibrium value of the lattice parame- 
ter of Ni in its bcc phase, introduce very little strain in the lattice. 
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TABLE 4: COORDINATION MATRICES FOR Ni, A1 AND Ti ATOMS IN TERNARY 
PHASES 


Heusler phase Other Ni 2 AlTi phase 



Nearest neighbors 

Next nearest 
neighbors 

Nearest neighbors 

Next nearest 
neighbors 

Atom 

Ni 

A1 

Ti 

Ni 

A1 

Ti 

Ni 

A1 

Ti 

Ni 

A1 

Ti 

Ni 

0 

4 

4 

6 

0 

0 

0 

4 

4 

6 

0 

0 

A1 

8 

0 

0 

0 

0 

6 

8 

0 

0 

0 

2 

4 

Ti 

8 

0 

0 

0 

6 

0 

8 

0 

0 

0 

4 

2 


Coordination of nearest-neighbor and next-nearest-neighbor bonds for Ni, A1 and Ti atoms in (a) the Heu- 
sler phase and (b) the alternate ordering pattern. 


A closer examination of the configurations in the Appendix shows that the preference for L2j 
ordering is also apparent in the low energy states in the nonselected configurations in Fig. 4. This 
is because any configuration that contains Ti atoms in A1 sites with only A1 next nearest neigh- 
bors introduces a significant change in energy. 




(b) Alternative phase A 2 BC 


Fig. 7: Schematic representation of (a) Heusler and (b) alternative ternary phases. 
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An additional advantage of this type of analysis is that a good deal of information on alterna- 


tive structures can also be obtained. This is an important issue when dealing with metastable 
structures that may result from a particular processing scheme. For instance, we could examine 
the energetics of an alternative type of ordering, which for example could involve a slight change 
in site occupancy in the Al- sublattice. Fig. 7 shows the Heusler unit cell as well as one corre- 
sponding to this alternative ordering scheme, which shares the first neighbor coordination with 
the Heusler phase. The only difference between these structures resides in the second neighbor 
coordination, with two Ti-Al bonds being replaced by pure Ti-Ti and Al-Al bonds, due to the dif- 
ferent distribution of Ti and Al atoms in the Al-sublattice. The equilibrium energy for each of the 
two cells is, within the precision of these calculations, identical. Therefore, the BFS strain energy 
contribution for each atom of a given species, in each cell, is the same. Moreover, within the BFS 
scheme, the only difference in energy between the two structures is due to the different type of 
next-nearest-neighbors bonds. The creation of Ti-Ti and Al-Al next-nearest neighbor bonds, pro- 
duces changes in chemical energy resulting in a noticeable change in total energy of formation 


TABLE 5 

BFS ENERGY CONTRIBUTIONS FOR ORDERED TERNARY PHASES 



Heusler phase 


Alternative phase 


Atom 

Strain 

energy 

(eV/atom) 

Chemical 

energy 

(eV/atom) 

Total BFS 
energy e 
(eV/atom) 

Strain 

energy 

(eV/atom) 

Chemical 

energy 

(eY/atom) 

Total BFS 
energy e 
(eV/atom) 

Ni 

0.108 

2.236 

1.929 

0.108 

2.236 

1.929 

Al 

0.503 

-1.128 

-1.240 

0.503 

-1.374 

-1.622 

Ti 

0.905 

-3.591 

-4.782 

0.905 

-2.680 

-3.338 


BFS strain (first column), chemical (second column) and total (third column) contributions for each of the 


two ordered stmctures shown in Fig. 7. 
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which favors the location of Ti and A1 atoms in the Heusler ordering scheme. Table 5 shows 
numerical results for the BFS strain, chemical and total energy terms for each type of atom in each 
type of structure. The total energy of formation for each of the two ordered structures is obtained 
by e = (2e Nj + e Aj + e Tj )/ 4 , where e x denotes the BFS energy contribution of an atom of spe- 
cies X. Consequently, the formation energy for the Heusler cell is -0.541 eV/atom whereas the 
alternating ordering scheme has a higher energy of -0.280 eV/atom. 

We then conclude that the changes introduced in chemical energy due to the creation of Ti-Ti 
and Al-Al next-nearest-neighbors bond makes this configuration energetically less favorable and 
therefore less likely to be found. However, depending on the processing conditions, it might not 
be improbable to find short range order similar to that of this alternate ordering scheme (Fig. 7.b) 
in addition to the expected Heusler ordering in rapidly solidified or quenched alloys. Moreover, in 
regions with low Ti concentrations, this phase competes with the Heusler phase in that they both 
share a great deal of common short-range order. 

The large difference in energy that can arise between various configurations with the same 
concentration is not reflected in the corresponding values of the lattice parameters (see Appen- 
dix). The lattice parameters of the different ternary alloys considered follow a nearly linear rela- 
tionship between the binary NiAl and NiTi B2 values as a function of Ti concentration (Fig. 8. a). 
There is, however, a small deviation above the ideal linear relationship, common to most of the 
configurations that display short- and long-range order. The overall trend, however, is clear. The 
lattice parameter of ternary Ni 5 oAl 5 o_ x Ti x alloys increases with increasing Ti content, as con- 
finned experimentally in Ref. 4. A more detailed display of these results is seen in Fig. 8.b, for Ti 
concentrations between 0 and 25 at. % Ti. 
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Table 6 lists the BFS predictions for the lattice parameter of the lowest energy configurations 
for each concentration. Two sets of values are shown: ' H easier ordering indicating the configu- 
ration where Ti atoms locate themselves exclusively in A1 sites (with only A1 next-nearest-neigh- 
bors) following Heusler ordering, and ' Solid Solution where the Ti atoms are randomly situated 
in the A1 sublattice. Note that there is no significant difference in lattice parameter for these two 
configurations. 

The results of the first column in Table 6 0 Ordered ’) can be easily adjusted to an expression 
of the form 

-2— = 1.00000 + 0.09347 x Ti (14) 

a NiAI 


TABLE 6 LATTICE PARAMETER OF 
Ni(Al,Ti) ALLOYS 


x Ti 

Heusler 

ordering 

Solid 

Solution 

2.78 

2.856 

2.856 

4.17 

2.859 

2.859 

5.55 

2.863 

2.863 

6.94 

2.867 

2.867 

8.33 

2.871 

2.870 

9.72 

2.875 

2.874 

11.11 

2.879 

2.877 


o 

Lattice parameter (in A) of the lowest energy NiAl-Ti configuration for several Ti concentrations (see 
Appendix for the lattice parameter values for the complete set of configurations used in this work). 
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plotted in Fig. 8.b, indicating the linear dependence of the lattice parameter of NiAl-Ti alloys with 
increasing Ti concentration, where a NjAi denotes the equilibrium value of the lattice parameter for 
stoichiometric NiAl alloys. 

The value of these results from this type of analysis is dependent on our ability to build a suf- 
ficiently complete catalogue of configurations that includes all the relevant atomic distributions. 
While the fit includes the lowest energy states that are part of the catalogue, no guarantee exists 
that these are actually the absolute lowest energy states for this system. 

The predicted values for the lattice parameter of NiAl-Ti alloys compare very well with exper- 
imental results by Kitabjian et al. [33]. In a recent study of atomic size effects of Ti in NiAl, they 
calculated the atomic size mismatch ali Ti for Ti atoms replacing Al atoms on the Al sublattice 


A Q r/ = 


Q Tj - Q 


Al 


Q 


Al 


( 15 ) 


where Q* n and n A/ are the atomic volumes of Al and Ti in NiAl. Using their calculated value Q A l= 
0.134 nm 3 , they obtain an approximate expression for Mi n , 


AQ r/ 


3 a p da 
2Q. A! dx T j 


( 16 ) 


where a 0 is the reference lattice parameter of NiAl at x Ni =0.5 and daldx Tj denotes the variation of 
the lattice parameter of NiAl-Ti alloys with Ti concentration. From measured values for a Ni- 
47.5Al-2.5Ti alloy, datdx Ti is estimated to be 0.0364 nm, so that the experimentally determined 
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value of a a Tj is 0.34, that is, the Ti atom is calculated to be 34 % larger than the A1 atom it 
replaces. 

The values of the lattice parameter as a function of Ti concentration predicted with BFS yield 
AQ Ti = 0.27 , also indicating a large size effect comparable with the experimental value. Consid- 
ering the constraints imposed in the theoretical calculations, where the physical properties of each 
element are computed from first-principles calculations [9], as well as the fact that the experimen- 
tal prediction is based on only one data point, the agreement between experiment and theory 
seems satisfactory. 

The BFS results are also in good agreement with the measured lattice parameters of NiAl-Ti 
alloys, including the Heusler alloy Ni 2 AlTi, reported by Oh-ishi et al. [34], Using their experi- 
mental values for NiAl-Ti and some alloys above the Heusler phase boundary (located at about 20 
at. % Ti), Kitabjian et al. [33] computed a size factor of 0.24 for Ti in NiAl, very similar to the 
value of 0.25 obtained from the corresponding BFS theoretical predictions (using Heusler phase 
values). 

The excellent agreement obtained by including Heusler phase values for the determination of 
size factors of alloying additions, encourages the use of this approach for those cases where the 
alloying addition has a very low solubility in the base alloy (i.e., Hf or Zr in NiAl) [33]. 

The BF Rule 

We can validate the theoretical value of these results in a general way by comparing them with 
alternative approximate algorithms for determining the lattice parameter. The lowest order 
approximation consists of using Vegard’s law [35], which establishes a direct linear average of the 
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lattice parameter of the constituents. However, this law is rarely valid even for the case of binary 
systems. Subsequent improvements on Vegard’s law highlight the role of compressibility in deter- 
mining intermediate values of the lattice parameter, in most cases introduced in an ad hoc fashion. 
Recently, a set of approximate expressions based on the BFS method was derived for predicting 
the lattice parameter of multicomponent systems which naturally introduces the compressibility 
of the constituents in its formulation. This expression, known as the BF rule has been proven in 
previous applications to multi-component systems [16], and provides an accurate estimate of the 
lattice parameter. 

The BF rule states that the lattice parameter of a multicomponent system (for elements of the 
same crystallographic species) is given by 

N 

S «1B; 

a(x I ,x 2 , ...) = (17) 

I «fi, 

i = 1 


where a i and B i are the equilibrium lattice parameter and bulk modulus, respectively, for the /th 
atomic species. Using the parameters listed in Table 1 yields the following expression for NiAl-Ti 
alloys: 


;V / ’ 


X AI’ 


X Ti ) 


1997.32* w + 794.43 jt a/ + 1249.44* rf 
685.80,v w + 248.88,r A/ + 388.87.vy. 


(18) 
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The results of Eq. 15 are also plotted in Fig. 8. Similar expressions (Eqs. (40)-(47) hi Ref. 16) 
can be obtained for the concentration dependence of the bulk modulus and cohesive energy per 
atom. For example, the bulk modulus for multicomponent alloys is given by the approximate 
expression 


B(x v x 2 , ...) = 



(19) 


Using the values of / from Table 1, together with its definition as a function of E c , a e and B 0 , 


I = 


j c 


127 iqa e B 0 


(20) 


we obtain the following approximate expression for the bulk modulus of NiAl-Ti alloys 


(685.80v Nj + 248.88a a/ + 388.87a n ) 2 
B ( x Ni’ x Al’ x Ti ) = 1997.32.1^ + 794.43.v /A/ + 1249.44.v y7 


(21) 


Also, simple mles can be obtained (Eqs. 45-47 in Ref. 16) regarding the deviation of such values 
to the usually assumed average values. Eq. 19 applies to any multicomponent system where all the 
constituents are parameterized in the same symmetry as the alloy under consideration, hi this 
work, Ni, A1 and Ti are described (Table 1) by their lattice parameter, cohesive energy and bulk 
modulus in the bcc phase. 
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Monte Carlo Simulations 


While the information provided by studying a large number of candidate configurations at 
absolute zero temperature provides valuable information on the energetics of the system at hand, 
it says very little in terms of the effect of temperature and processing on the micro structure of the 
alloy. In this sense, several numerical tools exist to investigate these issues. In this work, we con- 
centrate on the application of Monte Carlo procedures to perform simulations of Ni-Al-Ti alloys. 
When added to the results of the previous section, the simulation results should provide a more 
complete theoretical picture of the microstmctural behavior of these alloys. 

The Monte Carlo procedure employed in this work is a variant of that used to simulate the 
Ising system [36]. The computational cell used in all simulations consists of a large number of 
atoms (1024 in most cases shown in this paper) arranged on a body-centered-cubic lattice. Bound- 
ary effects are minimized through the use of periodic boundary conditions in all directions. The 
choice of the size of the cell results from a compromise between computational economy, conve- 
nience in visualizing the results, and the ability to distinguish between all the pertinent physical 
features of the system. Although the simulations shown correspond to a wide range of tempera- 
tures, for simplicity we ignore the lattice parameter dependence on temperature and use the zero 
temperature predictions discussed in the previous section, as we are mostly interested in ground 
state effects. 

For the (high) initial temperature in every simulation, the computational cell is a random 
solution, i.e. the Ni, A1 and Ti atoms are randomly assigned to each site, in direct proportion to 
their assigned composition (unless otherwise indicated). No vacancies were allowed in the major- 
ity of simulations reported in this work, though their incorporation does not present any added 
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degree of difficulty, as will be demonstrated shortly. Most of the simulations shown correspond to 
a process where a sequence of decreasing temperatures (the 'cascade 7 ) is chosen, where the sys- 
tem is allowed to equilibrate sequentially at each temperature. This simulates the 'slow cooling’ 
of the actual alloy, starting from a high temperature disordered solid solution. The equilibration 
procedure involves the random selection of a pair of atoms and their subsequent reversal in posi- 
tion. The reversal in chemical identity is accepted or rejected using the Metropolis criterion 
where, for a given temperature, an exchange is accepted if it lowers the energy or assigned a 
probability exp{-AE/kT), where k is Boltzmann constant and A E is the change in BFS energy 
between the configurations after and before the switch. The number of switches allowed is set as 
an input parameter and it determines, together with the difference in temperature between two 
successive steps in the simulation, the simulated cooling rate of the sample. A large number of 
switches can be considered equivalent to a longer stabilization time, which together with small 
changes hi temperatures simulate a slower cooling rate. In this calculation we ignored the depen- 
dence of the lattice constant with temperature, therefore temperature enters into the calculation 
only through the Metropolis criterion. After the system has achieved equilibrium (based on the 
total energy of the computational cell), various properties of the system are computed and aver- 
aged and then the simulation proceeds to the next temperature decrement. The properties calcu- 
lated include the average energy of the cell, the specific heat, and bond correlations. 

While these simulations do not attempt to mimic the detailed dynamics of the equilibration 
process, they do offer a qualitative view of the effects of rapid versus slow cooling of the system. 
The cooling rate ( that is, the size of the steps between the various temperatures considered in the 
cascade) is of critical importance in determining the final state of the system. Slow cooling results 
in a highly- ordered low-temperature state, while rapid cooling results in a more disordered mate- 
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Ni A1 Ti 

Ni 0.2256 0.3679 0.4065 
A1 0.7358 0.0771 0.1870 
Ti 0.8130 0.1870 0.0000 




Ni 


• A 1 

l\l 


Ti 




Ni A1 Ti 

Ni 0.0059 0.4966 0.4976 
A1 0.9932 0.0020 0.0049 
Ti 0.9951 0.0049 0.0000 



(a) (b) 

Fig. 9: Final structures of a Monte Carlo/Metropolis/BFS simulation on a 1024 atom cell of a Ni-25A1- 
25Ti alloy. Both (initially random) final states are obtained by lowering the temperature in different ways: 
Fig. 9. a shows the final geometry for a rapid cooling process, where the final temperature is reached by 
‘freezing’ the initial, high temperature, state. Fig. 9.b is obtained by slowly lowering the temperature in 
100 K intervals and allowing the cell to equilibrate, until the final temperature is reached (‘cascade’ pro- 
cess) at room temperature. The inset includes the coordination matrix for the final state of the simulation: 
the element ij in this matrix indicates the probability that an atom i has an atom of species j as a nearest- 
neighbor. The labels ij take the values 1, 2 and 3, corresponding to Ni, A1 and Ti. 
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rial often containing antiphase boundaries. As with actual processing, the temperature treatment 
of the sample is essential in determining the final state. It is to be expected that the slow cascade 
processes, used hi most of our calculations, will result hi highly ordered compounds, hi contrast, 
sudden cooling of the sample will result in regions of disorder often in the form of antiphase 
boundaries . It is also possible that additional phases - like the ones described in Fig. 7 and dis- 
cussed in the previous section - might also appear with a frequency proportional to the difference 
in energy of formation with respect to other structures (such information can be obtamed from the 
‘static’ analysis described in the previous section). The presence of these higher energy structures 
in the fast cooled sample can arise if quenching of the sample 'freezes’ specific domains within 
the cell, whose seed is already present in the initial disordered state. 

An example of the effect of cooling rate on structure is shown in Fig. 9 which includes the 
results of two separate simulations for an alloy with a bulk composition equivalent to Ni 2 AlTi, 
both starting from the same random state at high temperature. To facilitate the visual interpreta- 
tion of the results, the 1024-atom cubic cell is stretched along the <001> direction. The rapidly 
cooled cell, shown in Fig. 9. a, is characterized by order and disorder features. There is a clear 
trend towards separation of Ni-rich and Ti-Al alternating planes, as well as some indication of 
Heusler ordering in about half the sample. The lower half of the cell shows a noticeable degree of 
disorder and even the existence of some energetically unfavorable situations, mostly in regions of 
large Ti or A1 concentrations. Fig. 9.b, in contrast, shows the results of a cascade process that had 
finer steps between temperatures (slow cooled). Even though the final temperature was the same 
as in the previous case (room temperature), this cell displays almost perfect Heusler ordering. 
Most of the antisite defects present in this sample are eliminated after subsequent re-heating and 
slow re-cooling - a second cascade - which results in an almost ideal Heusler phase with a few 
antisite atoms. 
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Fig. 10: Final geometry for a cascade simulation (see text) 
on a Ni-49Al-lTi cell of 1024 atoms. An expanded view of 
the cell is also shown. 



Figure 10 displays similar results for the slow cooling of a Ni 50 (Al,Ti) 5 0 alloy with 1 at. % 
Ti. Starting with a random alloy at high temperature, the slowly cooled sample settles into a per- 
fect B2 NiAl ordering with Ti atoms in A1 sites. However, no trend toward precipitate formation 
is seen at this concentration, consistent with the configurational analysis described in the previ- 
ous section. 

Figure 1 1 shows results for a Ni 5 o(Al,Ti) 5 o alloy with 5 at. % Ti for two different random 
computational cells. The first three cells in Fig. 1 1 show the final results of three consecutive 
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1024-atom cell - Triple cascade 4608-atom cell 


Fig. 11: Results of three consecutive temperature cascades on a Ni-45Al-5Ti 1024-atom computational 
cell. The fourth column corresponds to a single cascade calculation on a larger (4608 atoms) cell. 
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temperature cascades on the same sample, each showing a higher degree of short range order rel- 
ative to the previous one, A multiple-cascade process consist of consecutive cascades where the 
final state of a given cascade becomes the initial state of the following one, except for the first cas- 
cade where the initial state is a random distribution of atoms. The last cell displays the results of a 
single temperature cascade for a large (4608 atoms) cell, showing essentially the same features 
observed in the smaller cell, thus ruling out inicrostructural effects based on cell size. Since this 
concentration is near the accepted solubility limit for Ti in NiAl, a larger cell allows for a closer 
examination of the final state and the possibility of precipitate formation, hi this case, the Ti atoms 
actually display three different types of behavior, as seen in the last column in Fig. 1 1, which cor- 
responds to the final temperature (T = 400 K): 1) the formation of Heusler precipitates, more 
clearly seen if the periodicity of the cell is taken into account, 2) Ti atoms in solid solution in the 
NiAl matrix with preference for the Al- sublattice, and 3) the distribution of Ti atoms following 
the alternative to Heusler ordering depicted in Fig. 7.b (atoms located in alternate sites along rows 
parallel to the edges of the sample) in regions often adjacent to Heusler precipitates. Also, a few 
Ni antisite atoms are seen. While it might prove to be a premature conclusion, it is noticeable that 
in the final state the antistructure Ni atoms seem to 'attract’ Ti atoms to their vicinity in the Al 
plane, creating a Ti-rich interphase with the NiAl matrix. This effect, while visible in the results 
of the simulation, could be also due to the limitation imposed by forcing the lattice parameter to 
remain constant during the simulation. This restriction could induce an unrealistic strain in the 
calculational cell that could eventually drive Ti atoms to the interphase in order to compensate for 
such effect. 

Similar structures are observed in a Ni-40Al-10Ti alloy which has undergone a double-cas- 
cade process (Fig. 12). In this type of process, an initially random cell is cooled to a given tamper- 
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Fig. 12: Final geometry for a double cascade process on a Ni-40Al-10Ti 1024-atom cell, (a) The first 
expanded cell results from being slowly cooled at equal temperature intervals from an arbitrary high tem- 
perature. (b) Then this cell is re-heated and slowly cooled again, reaching the final state shown in the sec- 
ond cell. 
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ature and then used as the initial state of a new cascade process. At the end of each cascade, the 
formation of Heusler precipitates (characterized by the chains of Ti atoms along <1 10> in the A1 
planes) is clearly observed. In spite of the re -heating and re-cooling of the sample, the Heusler 
precipitates reappear proving that the dominant effect of sufficient additions of Ti to Ni-rich NiAl 
is the formation of such precipitates. Some of the other features pointed out in Fig. 1 1 are also fea- 
tured in the 10 at. % Ti case: the solid solution of Ti atoms in the matrix, and the presence of an 
occassional antistructure Ni atom in the A1 sublattice, with an apparent clustering of Ti atoms in 
the vicinity. Also, some Ti atoms following the alternative ordering scheme previously discussed 
are clearly seen. 

Finally, we comment on the role of vacancies on the micro structure of the alloys studied. So 
far, all the examples shown ignore the presence of vacancies. If vacant sites are allowed in the cal- 
culation, it is found that none of the essential microstructural features discussed earlier are 
affected to any significant degree. However, the results show clear indication of vacancy cluster- 
ing in such a way that inner A1 surfaces are created. This is not surprising, given the large size of 
A1 atoms and the low surface energy of Al, both features conducive to the formation of A1 sur- 
faces. Moreover, due to the tendency of vacancies to coalesce, Ni antistructure atoms are then 
found in the vicinity of vacancy clusters. Some of the unfavorable Ni-Ni bonds thus created are 
compensated by the migration of some Ti atoms to Ni sites creating favorable Ti-Al bonds. 

EXPERIMENTAL ANALYSIS OF NiAl-Ti ALLOYS 

Three NiAl single crystal alloys (Ni-47Al-3Ti, Ni-45Al-5Ti and Ni-43Al-7Ti) were grown by 
a Bridgman technique at the University of Florida. The ingots were homogenized for 32 hrs. at 
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1644 K, aged for 6 hrs. at 1255 K, and slowly furnace cooled from the aging temperature. The 


purpose of this heat treatment was to produce a low temperature ‘equilibrium’ microstructure that 
would best correspond to the ground state conditions modeled under the BFS method. Samples 
for transmission electron microscopy (TEM) were prepared from 3 mm diameter cylinders elec- 
tro-discharge machined from the heat treated ingots. Slices sectioned from the cylinders were 



Fig. 13: (a) Bright-held TEM image of the microstructure of the Ni-47Al-3Ti alloy and (b) corresponding 
<110> SADP. The only features observed in the alloy are an occasional <100> dislocation marked ‘D’ on 
(a). 
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Fig. 14: (a) Bright-field and (b) dark-field TEM image of the Ni-45Al-5Ti alloy showing precipitation of 
fine Heusler particles, (c) the corresponding SADP and (d) an indexed, simulated pattern. 
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mechanically ground and electrochemically thinned in a twin-jet Tenupol-3 polisher. Microstruc- 


tural (bright-field/dark-field) and diffraction analysis were conducted in a Phillips 400T TEM 
equipped with a double tilt goniometer. 

Fig. 13. a shows a brightfield image of the Ni-47Al-3Ti alloy. Except for an occasional disloca- 
tion, the microstructure is very clean and featureless and shows no sign of any second phase pre- 




Fig. 15: (a) Dark- field TEM image of the Ni-43Al-7Ti alloy showing dense precipitation of rectangular- 
shaped Heusler precipitates and (b) the corresponding SADP . 
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cipitation. This is confirmed in the corresponding <110> zone-axis selected area diffraction 
pattern (SADP), presented in Fig. 13. b, which shows only the NiAl matrix spots and no extra dif- 
fraction features (spots or streaking) due to precipitation. 

In contrast, Fig. 14. a shows a bright-field image of the Ni-45Al-5Ti alloy after the same ther- 
mal treatment. Precipitation of a high density of extremely fine second-phase particles can be 
clearly seen, especially in the dark-field image shown in Fig. 14.b, where the precipitates appear 
bright on a black background. The corresponding <110> SADP in Fig. 14.c shows distinct extra 
spots which were indexed to a fee crystal structure with lattice parameter a 0 = 0.586 nm, corre- 
sponding to the Ni 2 AlTi (Heusler) phase. From the crystallographic information revealed from the 
diffraction pattern (schematically illustrated and labeled in Fig. 14.d), it can be seen that the 
Ni 2 AlTi phase nucleates with a cube-on-cube orientation relationship with the NiAl matrix, i.e. 
[110] N iAi//[110]Ni2AlTi an< i (001) N iAi//(001)Ni2AlTi- Dne to the fine size of the precipitates, 1-5 
nm, and the small lattice misfit between these two phases, on the order of 1.5 %, the precipitates 
are coherent with the matrix resulting in significant coherency strains around the particles, as seen 
in the bright-field image in Fig. 14. a. 

Figure 15 shows a dark-field TEM image of the Ni-43Al-7Ti alloy after the same heat treat- 
ment. Precipitation of a high density of well defined and coherent Heusler precipitate plates rang- 
ing in size between 10-50 nm can be clearly seen. While the precipitates are on average 10 times 
larger than those in the Ni-45Al-5Ti alloy, due to the small lattice misfit between the precipitate 
phase and the NiAl, most of the plates are still coherent with the matrix. Fig. 15.b shows the cor- 
responding SADP. 
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DISCUSSION 


In exploiting the computational simplicity of the BFS method in calculating the energetics of 
various alloy structures, it must be understood that the technique does not automatically provide 
the ground state configuration of a particular alloy but depends on the inclusion of that configura- 
tion in the catalogue of configurations selected for study. For reference in this and future studies, 
the present catalogue of configurations has been defined in the Appendix. Returning one last time 
to the issue of completeness of this set of configurations, Fig. 4 can be used as a clear example of 
how an occasional omission in the predetermined catalogue can be easily detected and corrected 
as the large size of the set provides enough statistical information to detect trends and patterns and 
therefore, infer what configurations might be missing. Had one of the energetically favorable con- 
figurations been left out of the original set, the omission would have been noticed as a 'disconti- 
nuity’ in the plot shown in Fig. 4.b. For example, had the Ni 2 AlTi Heusler phase (denoted with a 
circle in Fig. 4.b) not been included in the set, it would have been easily inferred by observing the 
series of states indicated with open circles, leading to the prediction of such a phase as the ground 
state for that concentration. In fact, the analytical BFS results for the 25 at. % Ti alloy not only 
suggest but confirm the L2 1 structure as the ground state for Ni 2 AlTi. Furthermore, the Monte 
Carlo-Metropolis simulations provide only the ground state configuration and serve as a comple- 
mentary technique and check on the analytical approach. 

The use of this survey method for studying the energetics of various alloy systems provides 
one with significant information on not only the ground state structure of a particular alloy but 
relates general trends in the energy of formation and lattice parameter to changes in concentration 
and atomic distribution. It also has the important additional advantage of being able to identify 
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metastable structures or configurations with energy close to that of the ground state, which may 
have a high probability of appearing in the alloy depending on the actual processing conditions. 
We note also that Monte Carlo simulations complement the use of our survey method to obtain 
ground state structures. 

One of the more significant results of this investigation is the successful application of the 
BFS method to the determination of the solubility level of a ternary addition to an ordered inter- 
metallic compound, in this case Ti in NiAl. The results of static calculations, shown in Fig. 4, sug- 
gested that the solubility of Ti in NiAl is close to 5 at. %. This was experimentally verified by 
examining the micro structure of three NiAl-Ti alloys, where it was found that Ni-47Al-3Ti was a 
complete solid- solution alloy with no second-phase precipitation, and the nucleation of a high 
density of Ni 2 AlTi precipitates appeared only in NiAl alloys containing 5 or more at. % Ti. 
Although the calculation of minimum energy ground state configurations via the BFS method was 
performed at 0 K and the experimental results were obtained at room temperature, the microstruc- 
ture of these high-melting point ordered alloys is not expected to be different at these low temper- 
atures. Moreover, the Monte Carlo simulations included in this work verify this point. Another 
useful result of the BFS calculations is the ability of the method to predict the atomic structure of 
the second phase particles and the lattice parameters of all the constituent phases. Both the static 
and Monte Carlo simulations were able to predict the correct atomic configuration of the resulting 
second phase particles, i.e. Ni 2 AlTi (Heusler phase), in NiAl-Ti alloys. Based on these calcula- 
tions, the lattice parameter of the Heusler phase was found to be 0.5828 nm, which is in close 
agreement with the average value of 0.5876 nm reported for the Ni 2 AlTi structure [36]. The lat- 
tice parameter of the corresponding solid solution NiAl-Ti alloy was calculated to be 0.2865 nm. 
Thus, based solely on the BFS analytical approach, the lattice mismatch between NiAl and 
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Ni 2 AlTi is found to be 1.7 %, which is very close to the approximately 1.5 % misfit obtained 
experimentally via TEM. A controlled mismatch between phases is one of the major design crite- 
ria in almost all high temperature alloys. Therefore, the ability to model not only the correct sec- 
ond phases but also the resulting lattice mismatch is a significant breakthrough in the 
computational design of high temperature alloys. 

An additional significant contribution of this analysis is the ease with which the site occu- 
pancy of a third element can be determined in a structure. Knowledge of site occupancy of an 
alloying element is an important and necessary piece of information in understanding the defect 
structure and its impact on mechanical properties [37]. In the case of NiAl-Ti alloys studied here, 
the BFS method was able to correctly predict the site occupancy of Ti in the NiAl lattice with 
additional information on the overall substitutional scheme, as seen in Fig. 6. The preference of a 
Ti atom to occupy an A1 site in Ni-rich NiAl is well established [4] and, in fact, all the empirical 
alloy development programs start with this basic assumption. 

The most common experimental techniques for determination of site occupancy are 
AFCHEMI (atom location by channeling enhanced microanalysis) [38] and APFIM (atom-probe 
field ion microscopy) [39], However, these are often tedious and very involved procedures and in 
the case of the fonner there are a number of complications which make the technique relatively 
inaccurate [40]. The ability and ease to analytically determine the site occupancy of alloying addi- 
tions, as proven by the BFS simulations shown in this work, is not only beneficial to an alloy 
developer, but is a boon to those trying to determine the mechanistic behavior of ordered alloys. 
Moreover, the methodology used in this work is easily applicable to more than one alloying ele- 
ment, which allows for a clear understanding of the interaction between the different alloying 
additions [41], as will be demonstrated in future articles. 
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Fig. 16: Summary of final states for Monte Carlo/Metropolis cascade processes for N^qA^q x Ti x alloys (x 
= 1, 5, 6, 7, 8, 10, 25) for 1024-atom cells. 


Ni-40Al-10Ti 


Ni-25Al-25Ti 
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While a great deal of structural information has been derived from the static calculations, they 
do not provide much information in regards to the effect of temperature and processing conditions 
on microstructure. Therefore, to complement the static calculations and further our understanding 
of these alloys, Monte Carlo procedures were used to perform simulations on the Ni-Al-Ti alloys. 
These results are summarized in Fig. 16. Figure 16 shows the 1024-atom computational cell for 
the final states of temperature cascades for a number of different concentrations, some of these 
have been discussed in detail previously. 

The results show that Heusler formation is apparent beyond the solubility limit for Ti. But as 
in any real system, there is also a statistical chance for the development of other structures at a 
rate based on differences in energy between them and the ground state. In this case, the alternative 
ordering scheme shown in Fig. 7.b and an occasional Ni- antistructure atom are evident. The 
abundance of this alternative type of ordering can be easily explained in terms of our earlier dis- 
cussion of nearest and next-nearest-neighbor coordination. With such small energy differences 
involved, when a Ti atom has the choice of occupying an A1 site in a Monte Carlo calculation, the 
difference between a Heusler site and a site belonging to the alternate ordering scheme makes the 
probability at high temperatures almost identical between the two sites. Therefore, once a Ti atom 
is 'trapped’ in an alternative ordering site, its likelihood to migrate further to form a Heusler 
arrangement with decreasing temperature diminishes. With increasing Ti concentration, the ener- 
getically favorable Heusler site becomes more prevalent and the formation of Heusler precipitates 
is clearly favored due to the larger number of energetically favorable bonds thus created. In addi- 
tion, the number of temperature decrements within the cascade - the cooling rate of the alloy - 
influences the density of these higher energy structures observed in the final micro structure. 
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CONCLUSIONS 


The BFS method has been successful in predicting the solubility limit, structure of the second 
phase particles, lattice mismatch between the alloy matrix and the precipitating phases, atomic 
size factors, and preferred site occupancy of the alloying additions. These results have proven that 
almost all the necessary parameters needed for a purely analytical alloy design approach are now 
within reach. The present results provide confidence in the BFS technique, the authenticity of the 
input parameters used (Tables 1 and 2) and the approach used for obtaining parameters by use of 
ab-initio methods, which removes the limitations imposed by the otherwise required experimental 
data base. The results also provide an energetic description of the detailed micro structure of the 
NiAl-Ti system. The real asset of the BFS method would be to model more complex systems such 
as ordered binary alloys containing two or more alloying additions, required for modelling alloys 
needed in practical applications, which is in progress. Work in this area will be presented in future 
papers. 
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APPENDIX 


The computational cell is defined in Fig. 3. a. It corresponds to a bcc lattice with 72 sites. The 
B2 NiAl alloy corresponds to the atomic distribution shown in Fig. 3. a, where Ni atoms are 
denoted by black disks (labeled 1, 2, 3, ...) and the A1 atoms are denoted by open circles (labeled 

13.14.15.. ..). A set of configurations is defined by changing the occupancy of these sites by 1) 
exchanging an A atom in site n with a B atom in site m (A <-> B) or by 2) substituting an atom B 
in site m with an atom A that originally was in site n (A -> B).Some of the configurations corre- 
spond to smaller versions of the 72 atom cell: those denoted with an asterisk (*) correspond to a 
cell where atoms 4n (n = 1,...,18) have been eliminated and those configurations denoted with a 
double asterisk (**) correspond to a fraction of the original cell where only atoms 
(1,2,5,6,9,10,17,18,25,26,29,30,37,38,41,42) are taken into account. Fig. 3.b shows two examples 
of the configurations included in the set listed below, corresponding to x Ti = 13.89 and x Ti = 

25.00. The Table in this Appendix lists, for each Ti concentration, the different configurations 
used and their resulting energy of formation and corresponding equilibrium lattice parameter 
(determined by minimizing the energy of formation). For example, one of the configurations for 
x Ti = 13.89 (shown in Fig. 3.b) is labeled as: [A1 -> Ti](37-48) + Ni(49,60) <-> Al(61,72), indicat- 
ing that the A1 atoms in sites 37 through 48 are replaced by Ti atoms, and that the Ni atoms origi- 
nally in sites 49 and 60, exchange places with the A1 atoms in sites 61 and 72. 

The other configuration shown in Fig. 3.b and corresponds to x Ti =25.00. The configuration is 
labeled as [A1 -> Ti] (14,16,17,19,22,24,37,39,42,44,45,47,62,64,65,67,70,72) and indicates that 
Ti atoms now occupy the listed sites in the Al-sublattice, so that the final configuration corre- 
sponds to the Heusler phase shown hi Fig. 7. a. 
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CATALOGUE OF ATOMIC CONFIGURATIONS 


x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

0.00 

B2 

-0.61311 

2.848 


Ni(31)<->A1(42) 

-0.51543 

2.857 


Ni(30,32)<->A1(4 1 ,43) 

-0.41837 

2.866 


Ni(30,32)<->A1(42,43) 

-0.41724 

2.867 


Ni(30,31)<->A1(42,43) 

-0.44122 

2.864 

1.39 

[Al<->Ti](42) 

-0.60954 

2.852 

2.78 

[Al<->Ti] (42,44) 

-0.60634 

2.856 


[Al<->Ti] (38,44) 

-0.60615 

2.856 


[Al<->Ti] (38,43) 

-0.60577 

2.856 


[Al<->Ti] (42,43) 

-0.60243 

2.856 

3.70 

[Al<->Ti] (39,66)* 

-0.60363 

2.858 


[Al<->Ti] (37,42)* 

-0.60350 

2.858 


[Al<->Ti] (42,43)* 

-0.59918 

2.859 

4.17 

[Al<->Ti]( 17,39,71) 

-0.60258 

2.859 


[Al<->Ti] ( 1 8,38,43) 

-0.60195 

2.859 


[Al<->Ti] (42-44) 

-0.59592 

2.860 


[Al<->Ti] (42,46,47) 

-0.59573 

2.860 


[Al<->Ti] (38,42,46) 

-0.59267 

2.860 


[Ni<->Ti] (50,54,58) + [Al<->Ti] (62,66,70) 

-0.41575 

2.873 

5.55 

[Al<->Ti]( 17,37,42)* 

-0.59861 

2.863 


[Al<->Ti] (37,42,47)* 

-0.59844 

2.863 


[Al<->Ti]( 15,18,38,43) 

-0.59816 

2.863 


[Al<->Ti](22,41,63,70) 

-0.59527 

2.863 


[ Al<->Ti] (41 ,42,47)* 

-0.59414 

2.864 


[Al<->Ti] (38,42,43)* 

-0.59048 

2.864 


NAS A/TM— 1999-208892 


58 










x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

5.55 

[Al<->Ti](38,39,41,42) 

-0.58915 

2.864 


[Al<->Ti]( 18, 19,42,43) 

-0.58637 

2.864 


[ Al<->Ti] (41 ,-43)* 

-0.58630 

2.864 


[Ni->Ti] (52,56) + [A1 -> Ni] (65-68) 

-0.35182 

2.882 

6.25 

[A1 -> Ti](13)** 

-0.59975 

2.865 

6.94 

[A1 -> Ti] (22, 4 1,64, 66, 72) 

-0.59194 

2.867 


[A1 -> Ti]( 15,18,38,43,46) 

-0.59098 

2.867 

7.41 

[A1 -> Ti] ( 14, 19,39,42)* 

-0.59383 

2.868 


[A1 -> Ti](14,42,45,67)* 

-0.59358 

2.868 


[A1 -> Ti](14,41,43,70)* 

-0.58945 

2.868 


[A1 -> Ti]( 15, 19,38,42)* 

-0.58562 

2.869 


[A1 -> Ti] (37,42,43,66)* 

-0.58560 

2.869 


[A1 -> Ti](17,42,65,67)* 

-0.58528 

2.869 


[Al -> Ti] (38, 4 1,43, 46)* 

-0.58466 

2.869 


LAI ->Ti](38,39,42,43)* 

-0.57847 

2.869 


[Al -> Ti](18,38,42,46)* 

-0.57803 

2.869 

8.33 

[Al ->Ti](15, 18,20,38,40,43) 

-0.59161 

2.870 


[Al -> Ti] (38, 40, 4 1,43 ,46, 48) 

-0.58482 

2.871 


[Al -> Ti] (20, 2 1,22, 46, 62, 68) 

-0.58235 

2.871 


[Al -> Ti] ( 17, 19,4 1 ,43,65,67) 

-0.57473 

2.872 


[Al ->Ti]( 17,41,42,43,44,65) 

-0.57093 

2.872 

9.26 

[Al ->Ti]( 14, 19,39,42,69)* 

-0.58892 

2.873 


[Al -> Ti] ( 14, 19,22,39,42)* 

-0.58465 

2.873 


[Al ->Ti] (37, 42, 45, 62, 70)* 

-0.58058 

2.874 


[Al ->Ti]( 14,22,37,39,42)* 

-0.58028 

2.874 


[Al ->Ti]( 14,37,39,42)* 

-0.57998 

2.874 
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x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

9.26 

[A1 -> Ti](38,39,42,43)* 

-0.57376 

2.874 


[A1 -> Ti] ( 17, 19,42,65,67)* 

-0.57312 

2.875 


[A1 -> Ti] (39, 42, 43, 45 ,46)* 

-0.57264 

2.875 


[A1 -> Ti] (38,41 ,42,43,46)* 

-0.56647 

2.875 

9.72 

[A1 -> Ti]( 13, 15, 18, 38, 41, 43, 46) 

-0.58443 

2.874 


[A1 -> Ti]( 17, 23, 39, 47, 6 1,66, 72 

-0.58159 

2.875 

11.11 

[A1 -> Ti](13, 15, 18,20,38,40,41,43) 

-0.58468 

2.877 


[A1 -> Ti] (38, 4 1,46, 6 1,66,69)* 

-0.57593 

2.878 


[A1 ->Ti](13, 18,38,41,43,66)* 

-0.57552 

2.878 


[A1 -> Ti](42,44,46,48,65,67,69,71)* 

-0.57302 

2.879 


[A1 -> Ti]((13,21,38,46,63,71)* 

-0.57235 

2.879 


[A1 ->Ti](17, 19,22,37,39,67)* 

-0.57197 

2.879 


[A1 -> Ti]( 13, 17, 19,42,45,66)* 

-0.57157 

2.879 


[Al ->Ti](14, 18,37,39,47,66)* 

-0.56802 

2.879 


LAI -> Ti]( 19, 20, 37, 39, 41, 44, 66, 67) 

-0.56679 

2.879 


[Al -> Ti]((15, 19,38,42,61,64,69,72) 

-0.56662 

2.879 


[Al ->Ti] (38, 4 1,42, 43 ,46)* 

-0.56176 

2.880 


[Al -> Ti] ( 13, 17,2 1 ,63,67,7 1)* 

-0.56116 

2.880 


[Al -> Ti] (37,38,42,43,45,47)* 

-0.56050 

2.880 


[Al -> Ti]( 13, 17, 18,2 1,42,66)* 

-0.55817 

2.880 


[Al -> Ti] (37, 39, 4 1,42, 43, 46)* 

-0.55792 

2.880 


[Al -> Ti](18, 19,22,23,42,43,46,47) 

-0.55335 

2.881 


[Al -> Ti] (4 1,42, 43, 44, 65, 66, 67, 68) 

-0.55220 

2.881 

12.50 

[Al -> Ti]( 13,38)** 

-0.58300 

2.881 


[Al -> Ti]( 13,42)** 

-0.58000 

2.881 


[Al -> Ti]( 13, 15,18,20,38,40,41,43,46) 

-0.57789 

2.881 
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x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

12.50 

[A1 -> Ti](13, 19, 22-24, 38, 43, 62,67) 

-0.56031 

2.883 


[A1 -> Ti] (13,14)** 

-0.55910 

2.883 


[A1 -> Ti] (38-40,42,44,46-48) 

-0.54184 

2.885 


[A1 -> Ti]( 14, 18,22,38,42,46,62,66,70) 

-0.53443 

2.886 

12.96 

[A1 -> Ti](38,41,43,46,63,66,69)* 

-0.57087 

2.883 


[A1 -> Ti]( 13, 15, 18,38,43,66,69)* 

-0.57081 

2.883 


[A1 -> Ti]( 18, 38, 4 1,46, 6 1,66,69)* 

-0.56697 

2.883 


[A1 -> Ti](13, 18,39,42,47,66,69)* 

-0.56350 

2.884 


[A1 -> Ti] ( 1 9,2 1 ,38,4 1 ,42,46,47)* 

-0.56090 

2.884 


[A1 -> Ti](18,37,39,42,45,47,66)* 

-0.55261 

2.885 


[A1 -> Ti](14, 18,22,37,38,46,47)* 

-0.55135 

2.885 


[A1 -> Ti]( 18,38,4 1-43,46)* 

-0.54795 

2.885 


[A1 -> Ti] (38,39,4 1-43,45,46)* 

-0.54360 

2.886 


[A1 -> Ti] (38,39,4 1-43,46,47)* 

-0.54099 

2.886 

13.89 

[A1 -> Ti](13, 15, 18, 20, 21, 38, 40, 41, 43,46) 

-0.57135 

2.885 


[A1 -> Ti](15, 18, 19, 38, 41, 42, 44, 47, 63,69) 

-0.55710 

2.887 

14.81 

[A1 -> Ti](l 4, 15, 17,22,37,43,62,71)* 

-0.56300 

2.888 


[A1 -> Ti]( 15, 18, 23, 38, 43, 6 1,63,66)* 

-0.55915 

2.888 


[A1 ->Ti](14, 17,19,22,37,39,42,45)* 

-0.55887 

2.888 


[A1 -> Ti]( 13, 17, 18, 22, 38, 4 1,42, 45)* 

-0.54781 

2.889 


[A1 -> Ti](18, 19,22,23,41,42,45,46)* 

-0.54140 

2.890 


[A1 -> Ti]( 13, 14, 17, 18, 37, 38, 4 1,42)* 

-0.53689 

2.890 


[A1 -> Ti]( 13, 14, 18, 22, 37, 38, 42, 46)* 

-0.53537 

2.891 

16.67 

LAI -> Ti] (37, 39, 42, 44, 45, 47 ,62, 64, 65,67 ,70, 72) 

-0.55922 

2.893 


[A1 -> Ti]( 15, 18, 23, 38, 41, 43, 6 1,66,71) 

-0.55407 

2.893 


[A1 -> Ti](14, 19,21,37,39,41,42,46,47)* 

-0.54840 

2.893 


[A1 -> Ti](38, 41, 43, 46, 6 1,63, 66, 69,71)* 

-0.54782 

2.893 
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x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

16.67 

[A1 -> Ti](14, 17, 22, 39, 42, 47, 61,65,69)* 

-0.54745 

2.893 


[A1 -> Ti]( 13, 23, 38, 39, 42, 45, 46,61, 71)* 

-0.54405 

2.894 


[A1 -> Ti]( 15, 19, 23, 37, 4 1,45 ,62, 66, 70)* 

-0.53747 

2.895 


[A1 -> Ti](39, 41, 42, 45, 46, 65, 66,69,70)* 

-0.53216 

2.895 


[A1 -> Ti](15, 19, 23, 37, 39, 41, 43, 45,47,61, 65, 69) 

-0.52729 

2.896 


[A1 -> Ti]( 18,37-39,42,45-47,66)* 

-0.52317 

2.896 


[A1 -> Ti](61-63, 65-67, 69-71)* 

-0.51218 

2.897 


[A1 -> Ti] (37-48) + Ni(49-60) <-> Al(61-72) 

-0.32156 

2.905 

18.52 

[A1 -> Ti](13, 18,21,38,39,41,62,67,69,70)* 

-0.54339 

2.898 


[A1 -> Ti](14, 17, 19, 42, 61, 63, 65,67, 69,71)* 

-0.53964 

2.899 


[A1 -> Ti](14, 17, 19, 22, 42, 61, 63, 66, 69,71)* 

-0.53938 

2.899 


[A1 -> Ti](13, 17,38,41-43,45,47,66,70)* 

-0.52873 

2.900 


[A1 -> Ti]( 18, 19, 38, 39, 42, 43, 63, 65, 66,71)* 

-0.52587 

2.900 

18.75 

[A1 ->Ti](13, 18,38)** 

-0.56750 

2.896 


[A1 -> Ti](13, 17,38)** 

-0.54137 

2.899 


[A1 -> Ti](13,14,18)** 

-0.52327 

2.901 

20.37 

[A1 -> Ti](14, 17, 19, 22, 37, 39, 42,45, 65, 67)* 

-0.52883 

2.904 


[A1 -> Ti](15, 21, 22, 37, 39, 41-43, 47,65-67)* 

-0.52345 

2.904 


[A1 -> Ti](15, 17, 18, 38, 41, 43, 46,62,65,67 ,70)* 

-0.52307 

2.904 


[A1 -> Ti](15, 18, 19, 23, 37, 38, 42, 46,61,67 ,70)* 

-0.52068 

2.904 


[A1 -> Ti](13, 15, 18, 37, 39, 41, 42, 45, 46, 61,66)* 

-0.51254 

2.905 

22.22 

[A1 -> Ti](14, 17, 19, 22, 37, 39, 42, 45, 47,62,67, 70)* 

-0.52146 

2.908 


[A1 -> Ti](13, 15,21,23,38,41,43,46,62,65,67,70)* 

-0.51472 

2.909 


LAI -> Ti](13, 18, 19, 21, 38, 39, 41, 42, 46,47 ,66, 67)* 

-0.51144 

2.909 


[A1 -> Ti](14, 17,18,23,38,39,41,43,45,47,63,67)* 

-0.50905 

2.909 


NAS A/TM— 1999-208892 


62 









x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

22.22 

[A1 -> Ti](13, 14, 19, 38, 41-43, 45, 62,66, 67,71)* 

-0.50631 

2.910 


[A1 -> Ti](13, 18, 37-39, 41, 43, 45-47,66,71)* 

-0.50456 

2.910 


[A1 -> Ti](13, 15, 17, 19, 21, 23,38, 42,46,62,66, 70)* 

-0.49771 

2.911 


[A1 -> Ti](13, 14, 17, 18, 21, 22,37,38,41,42,45, 46)* 

-0.48996 

2.911 

24.07 

[A1 -> Ti](13, 15, 18, 21, 38, 41, 43, 46,61,63, 66, 69,71)* 

-0.51303 

2.913 


[A1 -> Ti](14, 18, 19, 22, 37, 39, 41, 43,45,47, 62,66, 70)* 

-0.48968 

2.915 

25.00 

[A1 -> Ti](14, 17,37,42)** 

-0.55583 

2.911 


[A1 -> Ti](14, 16, 17, 19,22,24,37,39,42) + 
[A1 -> Ti](44, 45, 47, 62, 64, 65, 67,70,72) 

-0.52068 

2.914 


[A1 -> Ti](13, 14,37,42)** 

-0.50999 

2.915 


[A1 -> Ti](13, 14,41,42)** 

-0.50504 

2.916 


[A1 -> Ti]( 14,37,38,42)** 

-0.49365 

2.917 


[A1 -> Ti]( 14, 1 8,37,38)** 

-0.48963 

2.917 


[Al -> Ti]( 13, 14, 17,18)** 

-0.47424 

2.919 


[A1 -> Ni](14,17,37,42)**+[Ni -> Ti]( 1,6,26,29)** 

-0.36515 

2.913 


[ Al -> Ni]( 14,17,38,41)** + [Al -> Ti]( 13,1 8,37,42)** + 
[Ni -> Al](2,5,26,29)** 

-0.29494 

2.919 


[Al -> Ni](13, 18,37,42)** + [Al -> Ti]( 14,17,38,41)** + 
[Ni -> Al](l,2,5,6)** 

-0.29401 

2.919 


[Al -> Ni](13,18,37,42)** + [Ni -> Ti](25, 26,29, 30)** 

-0.23489 

2.925 


[Al -> Ni]( 17,18,41,42)** + [Al -> Ti]( 13,14,37,38)** + 
[Ni -> Al](5,6,29,30)** 

-0.23396 

2.925 


[Al ->Ti](13, 16, 17, 20, 21, 24, 37, 40, 41,44, 45, 48, 61, 65, 69) 
+ [Ni -> Ti] (27,3 1 ,35) + 

[Al -> Ni] (64,68,72)+ Ni(2, 6, 10, 50,51,54, 55, 58, 59) <-> 
Al( 14,18,22,38,42,46,62,66,70) 

-0.20455 

2.935 

25.92 

[Al -> Ti](13, 15, 18, 21, 23, 38, 41, 43, 46,61, 63,66, 69, 71)* 

-0.50165 

2.918 


[Al -> Ti](13,14,17,18,38)* 

-0.49371 

2.919 

31.25 

[Al -> Ti]( 13,18,37,38,41)** 

-0.48602 

2.930 


[Al -> Ti](13,14,17,18,38)** 

-0.44858 

2.934 
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x Ti 

Configuration 

AH 

(eV/atom) 

a 

(A) 

37.50 

[A1 -> Ti]( 14, 17, 37, 38, 4 1,42)** 

-0.42697 

2.947 


[A1 -> Ti]( 14, 17, 18, 37, 38, 41)** 

-0.42467 

2.948 


[A1 -> Ti] ( 17, 1 8,37,38,41 ,42)** 

-0.41253 

2.949 

43.75 

[A1 -> Ti] ( 14, 17, 1 8,37,38,4 1 ,42)** 

-0.38544 

2.963 

50.00 

[A1 -> Ti] ( 13, 14, 17, 1 8,37,38,4 1 ,42)** 

-0.35202 

2.976 
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